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Abstract 

The effects of pre- (i.e., gestation and during lactation) and post-weaning diet on the composition of faecal bacterial 
communities and adipose expression of key genes in the glucose and insulin pathways were investigated in the cat. Queens 
were maintained on a moderate protein:fat:carbohydrate kibbled ("Diet A"; 35:20:28% DM; n = 4) or high 
protein:fat:carbohydrate canned ("Diet B"; 45:37:2% DM; n = 3) diet throughout pregnancy and lactation. Offspring were 
weaned onto these diets in a nested design (n = 5 per treatment). Faecal samples were collected at wk 8 and 17 of age. 
DNA was isolated from faeces and bacterial 16S rRNA gene amplicons were analysed by pyrosequencing. RNA was extracted 
from blood (wk 18) and adipose tissue and ovarian/testicular tissues (wk 24) and gene expression levels determined using 
RT-qPCR. Differences (P<0.05) in composition of faecal bacteria were observed between pregnant queens fed Diet A or B. 
However, pre-weaning diet had little effect on faecal bacterial composition in weaned kittens. In contrast, post-weaning diet 
altered bacterial population profiles in the kittens. Increased (P<0.05) abundance of Firmicutes (77% vs 52% of total reads) 
and Actinobacteria (0.8% vs 0.2% of total reads), and decreased (P<0.05) abundance of Fusobacteria (1.6% vs 18.4% of total 
reads) were observed for kittens fed the Diet A compared to those fed Diet B post-weaning. Feeding Diet B pre-weaning 
increased (P<0.05) the expression levels of INRS, LEPT, PAI-1 and tended to increase GLUT1, while the expression levels of 
IRS-1 in blood increased in kittens fed Diet A pre-weaning. Post-weaning diet had no effect on expression levels of target 
genes. Correlations between the expression levels of genes involved in glucose and insulin pathways and faecal 
Bacteriodetes and Firmicutes phyla were identified. The reasons for why post-weaning diet affects microbial populations 
and not gene expression levels are of interest. 
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Introduction 

Mammalian health is driven by a complex interaction between 
the host genome, its diet and the composition and function of the 
microbiota that inhabit its intestine. The plasticity of the host 
genome is well established; it responds to changes in diet in part 
via epigenetic regulation of gene expression [1,2]. Pre-weaning 
diet (gestation and lactation) is known to affect the function of the 
genome and can have long term effects on the health of offspring 
[3,4] . While information is available on the effects of pre- and 
post-weaning diet on changes in gene expression [5,6], less 
information is available about these periods on the composition of 
intestinal microbiota. 

The stability of intestinal microbiota is of interest, due to the 
relationship between the diet and the metabolic health of the host 
[7] . The intestinal microbiota of adult cats can change in response 
to short-term diet [8,9], however it has been suggested that in the 



long term intestinal microbiota profiles appear to be stable [10]. 
Diet-induced-changes in microbiota have been linked to impaired 
intestinal function (e.g., increased intestinal permeability [11]) 
giving rise to disorders such as obesity [12], due in part to the 
ability of the intestinal microbiota to extract more energy from the 
diet and make it available to the host [13]. Therefore, many 
studies have focussed on the host at weaning, as the intestinal 
microbiota is still developing and is pliable [14]. 

More recently it has been recognised that the microbiota 
present at birth have long-term impacts on an individual's health 
and has a potential role in the risk of disease [15]. It is typically 
accepted that the intestinal tract of the newborn is sterile [16], 
however there is some evidence to suggest that there is 
translocation of microbiota from the mother to the offspring 
before birth [17,18]. Further, recent opinions suggest that 
maternal diet during gestation may influence the composition of 
the microbiota in the intestinal tract of the offspring at birth [19]. 
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Despite this, there is a lack of direct evidence to examine the 
effects of maternal diet on the composition of the offspring's 
microbiota, with only a few papers investigating how maternal 
nutrition impacts the microbiome of the offspring [20,21,22]. 

We hypothesised that pre-weaning diet (gestation and lactation) 
would alter the composition of faecal microbiota and the 
expression of genes in the glucose and insulin metabolic pathways 
offspring. In order to investigate the hypothesis, we selected diets 
that contained contrasting levels of protein, fat and carbohydrate, 
that have shown to promote large differences in faecal bacteria 
profiles of adult cats [9] . 

Materials and Methods 

Animal Study 

Ethics statement. The protocol for this study was approved 
by the Massey University Animal Ethics Committee (MUAEC 10/ 
108). All cats used are owned by Massey University and were 
housed at the Centre for Feline Nutrition (Massey University, 
Palmerston North, New Zealand) according to the Animal Welfare 
(Companion Cats) Code of Welfare (2007). 

Diets. Commercially available dry (kibbled; Diet A) and wet 
(canned; Diet B) diets were utilised in this study. Both diets were 
formulated to meet all nutrient requirements for growth, gestation 
and lactation according to the Association of American Feed 
Control Officials (AAFCO). Diets were analysed for moisture 
content using a convection oven at 105°C (AO AC 930.15, 925.10) 
and the ash residue using a furnace at 550°C (AOAC 942.05). 
Crude protein and crude fat were determined using the Leco total 
combustion method (AOAC 968.06) and acid hydrolysis/Mojon- 
nier extraction (AOAC 954.02), respectively. Gross energy (kj/ g) 
was determined using bomb calorimetry. Crude fibre was 
determined using the gravimetric method (AOAC 978.10) and 
Nitrogen Free Extractables (NFE) by difference (Table 1). 

Animals and housing. Eight queens (bred from 3 unrelated 
toms; mean age 5.3 ± 0.4 yr) were randomly allocated to one of 

Table 1. Macronutrient profile of a commercially available 
Association of American Feed Control Officials (AAFCO) feed 
tested maintenance diets fed to domestic short hair kittens 
(Felis catus). 



Component 


Diet A 1 


Diet B 2 


Dry Matter (DM; % as is) 


92.9 


31.7 


Crude Protein (% DM) 


35.3 


45.3 


Crude Fat (% DM) 


20.2 


37.6 


Ash (% DM) 


7.5 


7.2 


Crude Fibre (% DM) 


1.8 


1.5 


NFE 3 (% DM) 


28.2 


2.0 


Gross energy (kcal/100 g DM) 


507.0 


618.1 


Metabolisable energy 4 (kcal/100 g DM) 


393.8 


485.1 



1 1ngredient list of Diet A (from pack): Corn, chicken and chicken meal, chicken 
digest, maize gluten, chicken tallow, tuna meal, poultry and poultry meal, 
iodinised salt, vegetable oil. 

ingredient list of Diet B (from pack): Meat by-products and meat derived from 
chicken, lamb, beef, and mutton; gelling agent; minerals; vegetable oil, 
emulsifier; colouring; vitamins, chelating agents. 

3 Nitrogen free extractables calculated by difference (100 - crude protein - crude 
fat - crude fibre - ash). 

4 Determined using modified Atwater factors of: crude protein (3-5 kcal ME/g 
DM), crude fat (8-5 kcal ME/g DM), NFE (3-5 kcal ME/g DM). 
doi:1 0.1 371 /journal.pone.0080992.t001 



two diets (n = 4 per diet). The queens received either a Diet A or 
Diet B ad libitum throughout the study (Table 1). Each queen was 
mated with a single unrelated male fed either Diet A or Diet B. 
Once pregnancy was confirmed (approximately 7-8 wk of 
pregnancy), the queens were moved to individual, multi-layered 
digestibility cages, containing nesting boxes on an upper shelf [23]. 
One week before birth, a faecal sample was collected within 15 
min of defecation, snap frozen in liquid nitrogen and then stored at 
-85°C until analysis. 

In total, twenty kittens were born and were weighed within 
24 hours (h) of birth, and every other day for 2 weeks (wk), and 
then at weekly intervals until the end of the study. 

From wk 0 to 4, all of the kittens received milk from their dam 
exclusively. Half of each litter was randomly assigned (within sex) 
onto Diet A and half onto Diet B, forming four dietary treatment 
groups (A-A, A-B, B-B, B-A; Table 2). At 4 wks of age, the kittens 
were weaned onto solid food in a gradual manner, receiving both 
the appropriate solid food and milk until wk 8, when the kittens 
were fully weaned. During the weaning process (4 wks onwards), 
kittens were group housed according to post-weaning diet for 
meals and socialisation, returning to their dam only for access to 
milk. By 8wks of age, kittens were housed exclusively according to 
their post-weaning diet. Diet and water were available ad libitum 
daily to allow for normal growth (Figure 1). 

Pre-weaning diet has been defined as the combined effects of 
gestation and lactation and early weaning, and the post-weaning 
phase as the point where the kitten is no longer receiving milk 
from its mother. 

Faecal sampling and preparation of barcoded ampli- 
fication of bacterial 16S rRNA genes and amplicon 
pooling. At 8 and 17 wks of age, the kittens were housed 
individually for 24 h and a fresh faecal sample was collected within 
1 5 min of excretion, snap-frozen in liquid nitrogen and stored at — 
85°C. 

Nucleic acids were extracted from faeces (30 mg) with a 
previously described combined bead-beating and phenol/chloro- 
form protocol [24] . Fusion primers designed for GS FLX titanium 
pyrosequencing were used for microbiota analysis. Forward 
primers incorporating the Roche GS FLX adaptor A (5' -CCA 
TCT CAT CCC TGC GTG TCT CCG ACT CAG-3'), a unique 
12-base molecular barcode [25], a two-base linker sequence and 
bacteria specific sequence Ba9F (5 '-GAG TTT GAT CMT GGC 
TCA G-3') [26], and a reverse primer incorporating adapter B (5'- 
CCT ATC CCC TGT GTG CCT TGG CAG TCT CAG-3') 
and template specific sequence Ba515Rmodl (5'-CCG CGG 
CKG CTG GCA C-3') modified from Lane et al. [27] were used 

Table 2. Treatment groups for determining the effects of pre- 
weaning (in utero and during lactation) and post-weaning diet 
on intestinal microbiota in the domestic kitten (Felis catus; 
n = 5 per treatment). 







Queen diet 




Post-weaning diet 


Post-weaning diet 






Diet A 


Diet B 


Pre-weaning diet 


Diet A 


A-A 


B-A 






n = 3 females, n = 2 
males 


n = 3 females, n = 2 
males 




Diet B 


A-B 


B-B 






n = 3 females, n = 2 
males 


n = 4 females, n = 1 
male 





doi:1 0.1 371 /journal.pone.0080992.t002 
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2500 




Week 

Figure 1. The effects of Diet A or Diet B on growth rates of kittens. Data are reported as means ± SEM. 
doi:10.1371/journal.pone.0080992.g001 



to amplify the bacterial 16S rRNA gene. Each PCR reaction 
contained 40 ul of Taq PCR MasterMix (Qiagen, Hilden, 
Germany), 28 ul reverse primer (0.6 u,M) and 8 jal of barcoded 



forward primer (2 |iM). A 19 |il aliquot of the reaction mix was 
transferred into a sterile tube to serve as no-template negative 
control. The remaining 57 ul were spiked with 3 JJ.1 of template 



300 



200 



<2 150 



100 




1000 1500 
Sequences per sample 



Figure 2. The effects of Diet A or Diet B on faecal microbial diversity in pregnant queens (Felis catus). The rarefaction curves based on 
the Chaol diversity index (at 97% sequence identity cut-off) indicate that faecal bacterial communities of cats fed Diet A ( — ; n= 3 cats) were less 
diverse than those of cats fed Diet B (-; n = 4 cats). Data are reported as means ± SEM. 
doi:1 0.1 371 /journal.pone.0080992.g002 
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Table 3. The effect of diets on the bacterial phyla (% of total 
sequences) in faecal samples of pregnant female cats {Felis 



catus) fed Diet A (n 


= 3) or Diet B (n 


= 4). 








Phylum 


Diet A 


Diet B 










Mean 


Mean 


SEM 


P value 


FDR 


Actinobacteria 


0.21 


0.12 


0.07 


0.47 


1.00 


Bacteroidetes 


0.31 


27.39 


4.67 


0.06 


0.30 


Firmicutes 


99.26 


46.12 


4.34 


0.01 


0.04 


Fusobacteria 


0.14 


24.98 


2.18 


0.01 


0.05 


Proteobacteria 


0.03 


1.02 


0.27 


0.15 


0.74 



P value indicates ANOVA significance of rank transformed data and False 
discovery Rate {FDR) indicates multiple testing adjusted P value. 
doi:1 0.1 371 /joumal.pone.0080992.t003 



DNA at a concentration between 20 and 40 ng 1 . Amplifica- 
tion was performed using a Mastercycler proS (Eppendorf, 
Hamburg, Germany) with the following cycle parameters: initial 
denaturation at 95°C for 2 min, 30 cycles of denaturing (95°C, 20 
s), annealing (52°C, 20 s) and elongation (72°C, 1 min), and a final 
7-min extension at 72°C. Amplicons were purified using a High 
Pure PCR product purification kit (Roche Diagnostics, Mann- 
heim, Germany) and sent to Macrogen (Seoul, Republic of Korea) 
for pyrosequencing using the GS FLX titanium system (Roche 
Diagnostics, Mannheim, Germany). 

Tissue gene expression. At 18 wks of age a blood sample (0.5 
ml) was taken and RNA extracted using the Total RNA Minikit 
(Blood/cultured cells) from Geneaid® (DNature, Gisborne, NZ) 
according to manufacturer's instructions. At 24 wks of age, during 
neutering, samples of ovarian tissue (females only), testicular tissues 
(males only), and abdominal fat tissues (females only) were taken. 
The tissue samples were washed in ice-cold saline and then 
submersed in RNAlater® overnight, then stored at -20°C until 
analysis. RNA was extracted from abdominal fat and reproductive 
tissues using an RNeasy Plus Universal Mini Kit (QJAGEN Pty Ltd, 
VIC 3148, Australia) according to manufacturer's instructions. 

Expression of fatty acid synthase (FASN), Glucose transporter 1 
(GLUT1) and 4 (GLUT4), insulin receptor (INSR), leptin (LEPT), 
hormone sensitive lipase (LIPE), uncoupling protein 2 (UCP2), 
insulin receptor substrate 1 (IRS-1) and 2 (IRS-2), plasminogen 
activator inhibitor-1 (PAI-1), peroxisome proliferative activated 
receptor-y (PPARG) were determined in blood and tissue samples 
using RT-qPCR as described previously [6,14,28,29]. Briefly, total 
RNA was reverse transcribed using Applied Biosystems High 
Capacity RNA-to-cDNA kits (Applied Biosystems Inc., Foster 
City, CA) and rt-qPCR on a Rotor-Gene 6000 thermocycler 
(Qiagen). The expression of FASN, GLUT1, GLUT4, INSR, 
IRS1, IRS2, LEPT, LIPE, PAL-1, UCP2 and PPARG were 
normalised against expression of GAPDH. 

Statistical Analyses. Growth rates were determined utilizing 
repeated measures in the Mixed Model procedure of SAS using a 
two-way ANOVA with maternal and post-weaning diets as the 
main effects. The log transformation was utilized prior to analysis 
when necessary to maintain homogenous variance. No effect of 
gender was observed. Sequences were processed using the 
Quantitative Insights into Microbial Ecology (QIIME) version 
1.5 pipeline [30]. Sequences were denoised and chimeric 
sequences removed using the ampliconnoise.py wrapper. Se- 
quences passing quality checking were assigned to samples by their 
12 bp barcodes, which were then binned into operational 
taxonomic units (OTU) at a minimum pair-wise similarity of 
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Figure 3. The effects of pre- and post-weaning feeding of Diet A or Diet B on faecal microbial diversity in kittens (Felis catus). The 
rarefaction curves based on the Chaol diversity index (at 97% sequence identity cut-off) indicate that faecal bacterial communities of kittens fed Diet 
B (A-B and B-B; n = 10) were more diverse than those fed Diet A (A-A and B-A, n = 10). Pre-weaning diet did not have an effect on community diversity 
in the kitten. Data are reported as means ± SEM. 
doi:1 0.1 371 /journal.pone.0080992.g003 



97%. Taxonomy was assigned to each OTU using the Ribosomal 
Database Project classifier using an 80% confidence threshold 
[31]. Beta diversity between samples at the minimum sequence 
depth attained (1640 sequences) was compared using weighted and 
unweighted UniFrac phylogenetic distances. Two factor (pre- 
weaning diet X post-weaning diet) multivariate analyses of 
community phylogenetic diversity was performed using the 
ADONIS function in the R package Vegan. Gene expression 
data was analysed two-way ANOVA with maternal and post- 
weaning diets as the main effects in R 2.14.1. Rank transformed 
bacterial taxa proportions from phyla and genera that made up > 
0.05% of total bacteria in at least 5 samples were analysed by two- 
way ANOVA with maternal and post-weaning diets as the main 
effects. No effects of gender were observed. Results are reported as 
mean and associated pooled standard error of the mean (SEM) 
and were considered significant at P < 0.05 and a trend when 0.05 
< P < 0.10. Correlation heat maps between microbiota 
abundance at the phylum level and gene expression generated 
using the CCA package (Ignacio Gonzalez and Sebastien Dejean 
2009. CCA: Canonical correlation analysis. R package version 
1.2. http://CRAN.R-project.org/package = CCA). 

Results 

Influence of pre- and post-weaning diet on growth rates 
of kittens 

Prior to weaning (wks 0-4), bodyweight was higher (P < 0.05) in 
kittens exposed to Diet B. Post-weaning, kittens fed Diet B were 
heavier (P < 0.05) than kittens fed Diet A (Figure 1). 



Pyrosequencing of bacterial 16S rRNA gene barcoded ampli- 
cons resulted in a total of 120,520 sequences after denoising and 
chimera removal, with an average number of 2,564 sequences per 
sample (range = 1,640 to 4,813). The number of OTUs identified 
was 1968. Sequence length was on an average of 481 bp (range = 
201 to 530 bp). 

Influence of diet on faecal bacterial community 
composition of queens 

Diversity of the resident bacterial community in pregnant 
queens was increased by feeding Diet B, as shown by the Chaol 
diversity estimator (Figure 2). Diet significantly affected the 
proportion of phyla observed in the faeces of pregnant queens, 
with decreased Fusobacteria and Bacteroidetes and increased 
Firmicutes observed in queens fed with the Diet A (Table 3). In 
pregnant queens fed Diet A, over 99% of the bacteria observed 
belonged to the phylum Firmicutes. In contrast, faeces from 
queens on Diet B contained Firmicutes (46. 1 % of total reads), 
Bacteroidetes (25.7% of total reads) and Fusobacteria (25.4% of 
total reads) as the predominant phyla. In total, 29 bacterial genera 
were observed in the faeces of the pregnant queens. 

Influence of pre-weaning diet on faecal bacterial 
communities in kittens 

Pre-weaning diet did not affect the proportion of phyla in the 
kitten faeces (Table 4), however, the proportions of some taxa were 
altered (Table 5). Kittens from mothers fed the Diet A showed 
increased abundance of Solobacterium, Peptococcaceae and Clostridiaks 
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Figure 4. Principal Coordinate Analysis plot of weighted Unifrac phylogenetic distances showing the similarities between bacterial 
communities of queens fed Diet A or Diet B and their offspring fed Diet A (B-A or A-A) or Diet B (B-B or A-B) post-weaning. 

Percentage of variation captured by each component indicated on axes. 
doi:1 0.1 371 /journal.pone.0080992.g004 



and Clostridium (Table 5) and decreased abundance of Megamonas 
and Allisonella (Table 5). 

Influence of post-weaning diet format on faecal bacterial 
community composition of kittens 

There was no effect of sampling age on the microbial 
populations (8 or 17 wk, data not shown), therefore the results 
from these two periods were pooled. In kittens, the diversity of the 
bacterial population was affected by post-weaning diet (Figure 3). 
Principal Coordinate Analysis (PCoA) of weighted and unweighted 
Unifrac phylogenetic distances showed that the overall microbiota 
community structure was most similar between offspring fed with 
the same post-weaning diet, regardless of the mother's diet (Figure 
4 and 5). Multivariate analysis indicated that post weaning diet 
had a significant effect on community composition (P<0.001), 
while differences between pre-weaning diet tended towards 
significance (P = 0.07). No significant interaction between pre- 
and post-weaning diets were observed (P = 0.49). 

Sequence types representing five different phyla were observed 
in the faeces of kittens. Firmicutes and Bacteriodetes were the most 
abundant phyla observed in the kittens irrespective of diet, 
followed by Fusobacteria and Proteobacteria, with sequences 
belonging to the phylum Actinobacteria were the least abundant. 
Dietary format affected the proportions of phyla with Firmicutes 
(77% of total reads), Bacteroidetes (19% of total reads) and 
Proteobacteria (1.6% of total reads) being predominant in kittens 
exposed to Diet A. In contrast Firmicutes (54% of total reads), 
Bacteroidetes (24% of total reads) and Fusobacteria (19% of total 
reads) being predominant in kittens exposed to Diet B. Post- 
weaning diet significandy affected (P < 0.05) the proportion of 



phyla, with decreased Fusobacteria and increased Firmicutes and 
Actinobacteria for kittens fed Diet A (Table 4). 

In total, 19 bacterial families and 26 bacterial genera were 
identified in the faeces of kittens (Table 5). Post-weaning dietary 
format significantly affected 19 genera (Table 5). Predominant 
genera for the kittens fed Diet A were Streptococcus (18.9% of total 
reads), Prevotella (16.5% of total reads) and Megaspheara (12.7% of 
total reads). For kittens fed Diet B, the most dominant genera were 
unclassified genera belonging to the families Peptostreptococcaceae 
(26.7% of total reads) and Fusobacteriaceae (15.4% of total reads). 
Major shifts relating to post-weaning diet included decreased 
percentage of sequences for Bacteroides and Fusobacterium and 
increased percentage of sequences for Megasphaera, Streptococcus and 
Lactobacillus for kittens fed Diet A post-weaning (Table 5). 

Influence of pre- and post-weaning diet on tissue gene 
expression in kittens 

Maternal diet affected (P<0.05) the expression levels of LEPT, 
PAI-1 and INSR and tended (P<0.10) to affect the expression of 
GLUT1 (Table 6) in blood. In abdominal fat, expression levels of 
IRS-1 was affected by pre-weaning diet (P<0.05). There was no 
effect of pre-weaning diet in the expression levels of genes in 
reproductive tissues. There was no effect (P>0.05) of post-weaning 
diet on gene expression levels in the kittens fed either diet (Table 
6). Despite only limited numbers of genes showing differential 
expression due to dietary treatment, strong correlations were 
observed between the various members of the microbiota and gene 
expression in blood, fat and reproductive tissue (Figure 6). 
Abundance of Firmicutes was negatively correlated with expression 
of most genes examined in blood and fat, whereas abundance of 
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Figure 5. Principal Coordinate Analysis plot of unweighted Unifrac phylogenetic distances showing the similarities between 
bacterial communities of queens fed Diet A or Diet B and their offspring fed Diet A (B-A or A-A) or Diet B (B-B or A-B) post-weaning. 

Percentage of variation captured by each component indicated on axes. 
doi:1 0.1 371 /journal.pone.0080992.g005 



Bacteroidetes tended to be positively correlated with gene expression 
in blood and fat. In contrast, the correlation between Bacteroidetes 
and Firmicutes abundance and gene expression levels showed the 
opposite pattern in reproductive tissue, with the Bacteroidetes 
showing a negative correlation with gene expression, while the 
Firmicutes showed a positive correlation. The correlation profiles of 
the Actinobacteria abundance and gene expression levels were also 
most closely similar to that of the Firmicutes compared to other 
bacterial groups. 

Discussion 

We hypothesised that pre-weaning diet would significantly 
impact the composition of the faecal microbiota of the kitten. This 
study shows that despite high faecal microbial diversity in pregnant 
queens, the pre-weaning diet had no effect on faecal microbiota 
composition in offspring, but gene expression levels in blood and 
adipose tissues were altered by pre-weaning diet. This suggests that 
the maternal influences on gene expression are long lasting 
whereas its effects on microbial populations may be relatively 
short-term despite their relative stability over the lifetime of the 
individual [7,10]. 

Previous studies in rodents suggest that pre-weaning diet has an 
effect on the microbiota of the offspring, especially in early life 
[21,32]. The current study using next generation sequencing 
identified that pre-weaning diet had a small impact on the 
composition of the microbiota in the intestinal tract of the kittens, 
with post-weaning diets having a major role in influencing the 
composition of the faecal microbiota. This is evidenced by the 
large differences in the distribution of phyla between the pregnant 



queens and their offspring. In the current study, the pre-weaning 
period included both gestation and lactation and the transition of 
the offspring from milk to solid food. The pre-weaned offspring 
were therefore exposed to milk from their mother as well as their 
post-weaning diet. Bacterial species composition of the milk was 
determined using both culture-dependent and independent 
methodologies (Cookson, unpublished data) and suggested in- 
creased species diversity in milk from cats fed Diet A. Despite this 
we did not see many changes in intestinal microbial composition 
that were attributed to pre-weaning diet. In contrast, post-weaning 
diet had a major impact on the faecal microbiota of the analysed 
kittens. 

Kittens fed Diet A post weaning had increased abundance of 
Firmicutes and Actinobacteria and decreased abundance of 
Fusobacteria compared to kittens fed Diet B post-weaning, 
consistent with results observed in adult cats fed similar diets [9] 
and the queens in this study. The majority of studies investigating 
microbes in healthy mammals have identified Firmicutes and 
Bacteroidetes as the predominant bacterial populations present in 
the intestinal tracts. Our study corroborates these results. 
However, the predominant phyla detected in this study, Firmicutes 
and Bacteroidetes, differed to those observed previously [32] who 
reported Firmicutes and Actinobacteria as the two most abundant 
phyla. 

Actinobacteria showed a decreased abundance in kittens fed 
Diet B post-weaning. Included in the Actinobacteria are members 
of the genus Bifidobacterium and the family Coriobacteriaceae. Both of 
these taxa are associated with improved health in rodent models. 
For example, Bifidobacteria are thought to improve intestinal 
health [33], while Coriobacteriaceae are associated with decreased 
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blood glucose and cholesterol metabolism in rodents [34,35]. In 
the current study, Cottinsella was significandy reduced in the kittens 
fed Diet B. Decreased C. aerqfaciens has been associated with low 
carbohydrate diets in human subjects [36]. The overall proportion 
of Actinobacteria in the current study (averaging 0.6%) was much 
lower than that reported previously for kittens (28.5%) [32] and 
adult cats (16%) [9]) fed similar diets. Bifidobacteria were not 
detected in our study, contrary to a previous study investigating 
the faecal microbiome of kittens fed high-protein (52% DM) diets 
[32]. Results from our current study are in agreement with an 
earlier study from our laboratory using adult cats [9], but differ 
from previous studies [32], probably due to differences in 
analytical methodologies (primer design etc). 

Bacteroidetes are associated with both protein and GHO 
digestion [37]. In the current study there were no differences in the 
abundance of members of the Bacteroidetes associated with diet; 
this agrees with results in adult cats fed similar diets [9] . Within the 
Bacteroidetes, the abundance of Bacteroides was higher in kittens fed 
Diet B diet post-weaning. In the current study, the abundance of 
Bacteroides was also higher in kittens fed the Diet B post-weaning. 
This contrasts with results in adult cats which showed a decrease of 
Bacteroides associated with a similar diet [9]. Bacteroidete 
proportions were higher in the current study (c. 20%) compared 
to that observed previously in kittens fed moderate or high 
protein: carbohydrate dry diets (0.2% [32]). 

Kittens fed Diet B had lower abundance of Firmicutes 
compared to kittens fed Diet A. The abundance of Firmicutes in 
kittens fed Diet A post-weaning were similar to those reported in 
kittens of the same age fed a dry, moderate protein:CHO diet [32]. 
Within the phylum Firmicutes, large shifts in the abundances of 
Lactobacillaceae (Lactobacillus), Veillonellaceae (Megasphaera), Peptostrepto- 
coccaceae {Unclassified Peptostreptococcaceae), and Streptococcaceae (Strepto- 
coccus) were observed between diets. Lactobacillus populations were 
higher in kittens fed Diet A compared to kittens fed Diet B. 
Increased Lactobacillus on Diet A is in agreement with results in 
kittens fed a dry, moderate protein:CHO diets [32] and adult cats 
fed dry diets [9]. Similar to [32] who observed increased 
Megasphaera associated with a moderate proteimCHO diet, 
Megasphaera spp. were also more abundant in kittens fed Diet A 
in the current study. It is possible that the increased Megasphaera 
proportions were due to the increase in Lactobacillus spp. observed 
in the current study. In contrast, Peptostreptococcaceae proportions 
were higher in kittens on Diet B, which was in agreement with 
results in adult cats [9]. 

Fusobacteria are proteolytic bacteria [38] and are seen in low 
levels of in healthy cats (0.3-5% of sequences) [32,39,40]. They 
have been associated with high protein diets in kittens [32], adult 
cats [9] and dogs [41]. In the current study, the abundance of 
phylotypes belonging to the order Fusobacterium was higher in 
kittens fed Diet B versus Diet A post-weaning (15 vs. 0.9%, 
respectively). 

Overall, Proteobacteria proportions were not different between 
kittens fed Diet A or Diet B, with only Sutterella spp. increased in 
kittens fed Diet B post-weaning. Previous research has shown 
increased Proteobacteria associated with a dry, high proteinxar- 
bohydrate diet in kittens [32] and wet diets in adult cats [9]. The 
protein levels in the current study were similar (35% vs. 45% in 
Diets A and B, respectively) to that fed to adult cats (32-42%) [9]. 

In recent years there has been increasing focus on the effects of 
diet on the intestinal microbiota and the role they might play in 
the health of the host, including the domestic cat. However, much 
of the research carried out in the domestic cat has been conducted 
using research-formulated diets, which are generally dry in nature. 
In the current study, the macronutrient contents of Diets A and B 
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Figure 6. Correlation heatmaps of faecal microbiota and gene expression in (A) blood, (B) adipose fat tissue, and (C) reproductive 
tissue. Cell colour indicates strength and direction of correlation, with blue representing maximal negative correlation and red representing maximal 
positive correlation. 
doi:1 0.1 371 /journal.pone.0080992.g006 



differed largely in fat (37.6% vs. 20.2% respectively) and 
carbohydrate (2.0% vs. 28.2% respectively) with protein being 
slightly higher in Diet B (45.3% vs. 35.3% respectively). Because of 
these differences we are unable to make any conclusions about the 
type(s) of nutrient(s) that might have caused the shifts in bacterial 
abundances observed, however, macronutrient digestibility (i.e., 
protein, fat) was decreased in kittens on Diet A (data not shown). 
However, in rodent models Proteobacteria and Firmicutes 
increased in association with high fat diets [42]. The diets in the 
current study also differed in their protein source. In rodent 
models, beef and milk (standard chow) proteins affected the 
phylogenetic composition and diversity of intestinal bacterial 
communities [43,44], with increased diversity associated with the 
beef-based diets. Increased bacterial diversity has been observed in 
adult cats fed a meat-based diet (Diet B) [9]. It is therefore 
plausible that protein source may also affect bacterial community 
composition in the cat. The different effects of high protein diets 
on Fusobacteria and Lactobacillus spp. in different hosts may therefore 
reflect evolutionary differences in diet preference and is of interest 
for future research. 

Due in part to the rise of obesity in the domestic cats, there is 
increasing literature on the effects of obesity [45] and insulin 
sensitivity [28,46,47] on gene expression in the domestic cat. 
Obesity in the domestic cat has been linked to decreased IRS-1 
and IRS-2 expression levels in liver and skeletal muscle [28] . IRS- 
1 is involved in insulin resistance and is involved in the regulation 
of glucose uptake and conversion into fat cells [47]. While obesity 
is also linked with decreased GLUT4 protein expression levels [45] 
lipid-infused cats show increased GLUT4 mRNA expression levels 
[46] . While bodyweight gain rather than diet is thought to increase 
the occurrence of pre-diabetic state in the domestic cat [48], the 
effects of maternal diet on gene expression in the domestic cat are 
of interest due to the long-term effects that these changes may have 
in relation to the health outcomes of the offspring [6]. In the 
current study, offspring fed Diet B pre-weaning increased the 
expression levels of many genes including those that are involved 
in the glucose (GLUT1, INSR) and lipid (LEPT, PAI-1) metabolic 
pathways and decreased IRS-1 expression levels in abdominal fat. 
Previous studies in kittens have shown that high protein, pre- 
weaning diets increases the expression levels of key genes in the 



insulin (e.g., IR), glucose (e.g., uncoupling protein 2) and lipid (e.g., 
LEPT) metabolic pathways [6]. 

There is increasing interest in the interaction between the host 
and its intestinal microbiota profiles and how they interact. In this 
study, Firmicutes and Bacteroidetes had opposite correlations with 
expression of genes in the glucose, insulin and lipid metabolic 
pathways. For example, in blood, Firmicutes are negatively 
correlated to glucose and insulin genes, whereas Bacteriodetes 
are positively correlated. However, in reproductive tissues 
Firmicutes are positively correlated to glucose and insulin genes 
and Bacteriodetes are negatively correlated. It is established in 
rodent and human that the ratio of Firmicutes:Bacteriodetes 
increases in obese humans and rodents [12,13,49,50]. This ratio 
has been associated in perturbations in energy homeostasis and the 
development of insulin resistance (see [51,52] and may explain the 
observations of the current study. In abdominal fat tissue, 
fusobacteria abundance was strongly correlated with LEPT, 
FASN, PPARG, PAI-1 and IRS-2 gene expression levels while it 
is was negatively correlated with LIPE. This suggests that 
Fusobacteria could also be a potential biomarker of perturbations 
in energy homeostasis and warrants further investigation. 
Furthermore, the Firmicutes and Actinobacteria had similar correla- 
tions with gene expression, suggesting overlapping or synergistic 
functions of these two taxa in the microbial community. 

In conclusion this study has identified that while pre-weaning 
diet plays a role in the alteration of gene expression levels it has 
little impact on the faecal bacterial profiles of the offspring. In 
contrast, while post-weaning diet format has no effect on gene 
expression levels in offspring, it plays a critical role in faecal 
bacterial composition. Why this may be the case is of interest for 
future research. Relationships exist between intestinal microbiota 
profiles and genes involved in the glucose and insulin pathways 
and these may play a role in long-term health. 
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